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In this paper we compute the change in the energy of a uniformly 
magnetized uniaxial platelet produced by the introduction of a cylin- 
drical domain. Differentiation, of the energy expression yields the 
translational force produced by gradients in plate thickness, material 
composition, or temperature. The force expressions provide a means for 
estimating the effect of gradients in these parameters on the margins of 
domain devices. Equating the sum of the gradient produced forces to 
the drag force yields a general domain velocity expression. The various 
results are presented in both graphical and tabular form. 

I. INTRODUCTION 

Magnetic memory and logic devices employing cylindrical domains 
in uniaxial platelets have recently received considerable attention. 1,2 
The theory of the static stability of these domains 3 and its application 
to cylindrical domain devices*-' have been discussed in previous papers. 
This paper is concerned with the translational forces acting on the 
domains and their effect on device performance. 

* Present address Department of Electrical Engineering, McMaster University, 
Hamilton, Ontario. 
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In most device applications, cylindrical domains are propagated by 
gradients in the applied field. 1 - 2 Cylindrical domains may, however, be 
propagated by gradients in any of the independent parameters which 
determine the total domain energy. These parameters are: the applied 
field, H; the plate thickness, h; the saturation magnetization, M g ; 
and the wall energy density, tr w . The domain radius, r. n , is not an 
independent parameter for domains in equilibrium but is determined 
once the other parameters are specified. The gradients in a w and M s 
may be produced by composition gradients, strain gradients or tem- 
perature gradients. Composition or thickness gradients may be used to 
provide forces which are functions of position only, while temperature 
or strain gradients may be used to provide time variable forces. 

The translational force is obtained by differentiation of the total 
domain energy expression with respect to position, under the assump- 
tion that the gradients in the domain parameters which produce the 
translational force are sufficiently small that the domain remains 
circular and stable ; consequently the energy expression remains valid. 
Since this method of computing the force is independent of the de- 
tailed stress pattern which produces the domain motion, no estimate 
of the shape distorting tendency of the various parameters is obtained. 
Equating the sum of the translational forces to the drag force yields 
the general velocity equation, and comparing the magnitudes of the 
various forces yields their effects on device operation. 

II. DOMAIN ENERGY 

The energy change produced by the introduction of a single isolated 
circular 180° domain into an infinite plate of uniaxial magnetic 
material which is otherwise uniformly magnetized along the average 
plate normal (the z axis) is now calculated. Such a domain configura- 
tion is shown in Fig. 1. The assumptions and notation of Refs. 3, 4 and 
5 are maintained except that the domain parameters h, H, M 8 , a w , and 
r are allowed to be functions of position on the plate. In particular, 
the following is assumed in the model: the wall has negligible width, 
the wall is everywhere parallel to the z axis, and the wall energy 
density is independent of wall orientation or curvature. The values of 
all parameters are assumed to vary sufficiently slowly that they may 
be represented by their 2-averaged values at the center of the circular 
domain. Additionally, the applied field is represented by its z com- 
ponent, H, since under the assumptions stated above only this com- 
ponent interacts with the magnetization. 
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Fig. 1 — The cylindrical domain configuration and coordinate system. 

The total change in the energy of the system due to the presence of 
a cylindrical domain is 

E T = 2vrohff„ + 2M,Hivrlh - 2irh 3 (2TrM])I(2r /h). (1) 

In this expression the first term, the wall energy, is the product of 
the wall energy density and the wall area ; the second term, the applied 
field interaction energy, is the product of the magnetization change, 
2M 8 . the applied field and the domain volume; and the third term, 
the internal magnetostatic energy, is the negative of the integral of 
the generalized radial magnetostatic force of Ref. 3. The internal mag- 
netostatic energy function, I{2r lt /h) is therefore defined as 



I(2r n /h) = / F(x) dx, 



(2) 



where F{x) is the force function defined by equations 33 and 138 of 
Ref. 3. The lower limit of the integral (2) is chosen so that when the 
plate is uniformly magnetized, /•<> = 0, the domain energy expression 
(1) is zero. Various closed form and power series representations of 
I(d/h) are given in the appendix of this paper. This function, which 
is plotted in Fig. 2 as a function of the diameter-to-thickness ratio, 
d/h, and is tabulated in Table I, has the asymptotic forms 



• 
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Fig. 2 — The internal magnetostatic energy function, /, as a function of the 
diameter-to-thickness ratio. 



*©«# -£+&©■ + 



b+m^i !«> « 



and 



'(! 



l(d 



2_ /d\ 3 J_ /dV 
3tt W + 16 U, 



(«- 



(3b) 



The physical origin of the asymptotic behavior of the / function was 
discussed in some detail in Section IV of Ref. 3. 
The domain energy per unit wall length is 



Er/rd = (2rM 2 ,)h 2 [l/h - J(d/h)], 



(4) 



where I = a w /^T?M 2 , is the characteristic material length, where the 
applied field has been eliminated using the equilibrium condition 3- ' 



and where 



1 _i_ d H - F\-\ o 
h^hlirM. \h) ' 



(5) 



(6a) 
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Table I — Magnetostatic Energy, Force and Stability Functions 
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The function J has the asymptotic forms 
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The function J(d/h), the normalized total magnetostatic energy per 
unit wall length, is plotted in Fig. 3 together with the force function, 
F(d/h) and the stability functions S (d/h) and S^id/h) which have 
the asymptotic forms 
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Numerical values of all these functions are given in Table I. 

Since from the figure and the asymptotic forms of the functions, 
(6), (8) and (9), J(d/h) lies roughly midway (with respect to diam- 
eter) between S (d/h) and &>{d/h) and since the condition for domain 
stability is 3 " 5 S {d/h) > l/h > S 2 (d/h), then a platelet of arbitrary 
thickness may always be biased such that the introduction of a domain 



CYLINDRICAL MAGNETIC DOMAINS 



717 



1.25- 



O >- 

<r t 

M 

o < 



0.50 



0.25 





/ 
















/ 
















/ 












^^ 




/ 
/ 




F^. 










/ 








S 2- 








/ 
















/ , 










J_ 






















/, 












"^r 




^ 

















4 

d/h 



Fig. 3 — The magnetostatic rncrgy per unit wall length function, J, the mag- 
netostatic force function, F, and the magnetostatic stability functions, So and & , 
as functions of the diameter-to-thickness ratio. 

produces either a positive, or negative, or zero change in the total 
energy. The force function, Fid/h), was included in the figure so that 
the bias fields yielding these conditions may be determined. The values 
of the diameters and bias fields of zero energy determined using Fig. 3 
or Table I together with (4) and (5) are found to agree with those 
obtained previously. In the case of a domain having the preferred 
dimensions, 4,5 l/h = 0.2500 and d/h = 2.000 (corresponding to an ap- 
plied field of H/4*M B = 0.279), then J (2.000) = 0.2422. Since J is 
nearly equal to l/h the total energy is nearly zero. Under the preferred 
conditions and in a platelet in which 4wM 8 = 100 Gauss and d = 10 
microns, the absolute value of the wall energy and the total mag- 
netostatic energy change are each approximately 0.2 times the rest 
energy of the electron. [Note that the terms in (4) may not be identi- 
fied as wall energy and magnetostatic contributions because the 
equilibrium condition was used to eliminate the applied field.] 



III. THE TRANSLATIONAL FORCE 

The translational force is given by 
F = -VE T , 
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In this expression h, H, M, and cr„ are considered to be independent 
variables (functions of position on the platelet) and r Q is a dependent 
function of these variables determined by the equilibrium condition (5) . 
Since the fundamental' equilibrium condition is (dE T /dr ) h , H .«„,M. = 0, 
the last term in equation (10) may be dropped. Evaluating the remain- 
ing terms using equations (1) and (2) yields 

F = - [27rr <r„, + 2M ,Urr\ - ^h\2TrM 2 .)I{2r Q /h) 
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Eliminating the applied field using the equilibrium condition (5) and 
rearranging yields 

F = -k dh\2irM\)\-[l/h - F h (d/h)]Vh/h + [l/h - F H (d/h)]VH/H 

+ [l/h + F M (d/h)]VM./M. - 2[l/h]V<T u ,/a w ), (12) 



where 



4) = #(!'- *>■(") = ■■'■'(? 



2d 



3_ 

32 



©* + [-!©> 



■k \h 



Z(d 

8\h 



!«•• 



and 



fVm-41,)/ 



M 



F H (d/h) = F(d//i), 



=i-l+k($+b+m 



In 



= d _ 2_ (W J_ (d\* 

h 3tt W + 64 \hJ ' 



i«- 



(13a) 



5«1, (13b) 



(13c) 
(14) 

(15a) 



J«C1, (15b) 
(15c) 



CYLINDRICAL MAGNETIC DOMAINS 



719 



Figure 4 shows plots of F h , F M and F = F H as functions of the norm- 
alized domain diameter. Numerical values of these functions are in- 
cluded in Table I. From the figure and the asymptotic forms of (7), 
(8), (13) and (14), it is seen that the functions are positive and that 
both F h and F H are greater than S (d/h). From these properties and 
the stability requirement S {d/h) > l/h > S 2 (d/h), it is seen that for 
any stable domain the thickness gradient and magnetization gradient 
contributions to the force are in the direction of the gradient while 
the field gradient and wall energy gradient contributions to the force 
are in the direction opposite to the gradient. 

The absolute value of each of the terms in equation (12) is small but 
measureable. If # and Ho are defined as the collapse and elliptical run- 
out fields respectively, then in the case of a gradient in the applied field 
where the stability limits arc known roughly (see Ref . 5) , the force pro- 
duced on a domain for which d/h = 2, h/l = 4, V// = (H„ - H 2 )/d, 
d = 100 microns and 4*M a = 100 Gauss is approximately 6 X 10" 3 
dynes. While such absolute force measurements could possibly be 
carried out, a more relevant experiment for device applications is the 
balancing of field gradients against gradients in the other quantities. 
At the present time, no such measurements have been completed. How- 
ever, the directions of the applied field gradient force, the wall energy 
gradient force, and the thickness gradient force have been verified. The 
sign of the H gradient force is verified in everyday device operation. 
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Fig. 4 — The transverse magnetostatic force functions, F u , Fk , and Fh (F = Fh) 
and the magnetostatic stability functions So and S 2 , as functions of the diameter- 
to-thickness ratio. 



720 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1971 

The sign of a w gradient force was checked using a temperature gradient 
in a sample of Smo.ssTbo.isFeOs . According to the data of F. C. 
Rossol, 7 at room temperature the wall energy temperature coefficient 
is approximately +1.5%/C° and the magnetization is only a slightly 
decreasing function of temperature. The <r w gradient force thus domi- 
nates in this material, and we observed that the domains move 
away from heated regions towards cooler regions as predicted. The 
sign of the M a gradient force was verified using a temperature 
gradient in a sample of Gdo^Tbo^FesO^ garnet. In this material M a 
increases at 3%/C° at room temperature while the wall energy is 
approximately constant. The M a gradient force thus dominates, and we 
observed that the domains moved towards heated regions as predicted. 
The force direction measurements were completed in tapered platelets 
of orthoferrite in which we observed that the domains move towards 
the thick end of a platelet when restraints were removed. 

It should be noted that in carrying out these experiments it is 
important to obtain low coercivity materials. This is especially true 
in measuring the thickness gradient force when Vh/h is small. 

The drag force acting on a domain propagating with uniform 
velocity v d from equation 58 of Ref . 5 is given by 

F d = -= dh M.\* H c + -?- |v d |] T^T , (16) 

where H c is the wall motion coercivity and p. w is the wall motion 
mobility. Equating to zero the sum of the gradient force (12) and the 
drag force (16) yields the velocity equation for an otherwise freely 
propagating domain. In order to avoid the vector sum in this equation, 
it is convenient to assume that all the gradients lie in the same direc- 
tion and are positive. It is also convenient to assume that the gradients 
are uniform so that their magnitudes may be expressed in terms of 
the maximum parameter difference across the domain divided by the 
domain diameter, ^X = &X/d. Under these assumptions the velocity 
equation is 

8 H c 2 \v d \ _ r Ah r AH „ AM, „ A^ , 

7r47rM. + hAtM. ~ th h Lh H + Cm M. L ° a w ' U7a) 

where 

h 
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and 



are called transverse force coefficients and are all positive. 

For a domain having the diameter d = Sim a plate of thickness 41, 
the velocity equation becomes 

8 H e 2 |vj 



= 0.238 f - 0.279 ^ + 0.772 ^ - 0.250 ^ (18) 

If the magnitude of the sum of the gradients is not sufficiently large, 
no motion takes place. For example, when a domain having dimensions 
such that equation (18) applies is subjected to a field gradient, the 
magnitude of the field difference across the domain must satisfy the 
condition AH/H > 9.1 X H r /4irM s before any motion can take place. 

The transverse force coefficients are plotted as functions of the 
normalized thickness h/l in Fig. 5 for the bias condition d = (dnd^)* 
where d n and d 2 are the collapse and elliptical runout diameters 
respectively. For small thicknesses the asymptotic forms of C h , Cm , 
and C a are proportional to (l/h) exp (—n-l/h) and the asymptote of Cn 
is proportional to exp (—n-l/h). For large thicknesses C H and Cm 
approach unity, and C a and C,, approach the asymptote (8/3^)* {l/h) i . 

Some caution must be exercised in interpreting equation (17) and 
Fig. 5. First, the stability of moving domains has only been investi- 
gated for the case of gradients in the applied field and then only 
incompletely (see Ref. 5). Another problem is that drive gradients 
which are applied from the surface and which must obey Laplace's 
equation, such as field gradients and temperature gradients, decrease 
exponentially into the platelet. (This does not apply to volume heating 
such as laser heating.) For any given value of these gradients at the 
surface the maximum z averaged value of these gradients which may 
be applied thus decreases according to the inverse first power of the 
plate thickness in thick plates. This consideration shows that the use 
of platelets having a thickness no greater than the preferred thickness 
of 41 is thus more strongly preferred for achieving a high domain 
velocity than is indicated by inspection of Fig. 5 alone. 
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Fig. 5 — The transverse force coefficients d, , Cn , Cu , and C a for the bias condi- 
tion d = (dodv)K as functions of the normalized thickness, h/l. 



IV. SIGNIFICANCE FOR DEVICE APPLICATIONS 

Equation (17) and Fig. 5, when properly interpreted, show again 
that there is a preferred plate thickness for achieving high bit rates in 
devices and that this thickness is again 41. Equation (18) shows that, 
for plates of this thickness, wall energy, magnetization and applied 
field gradients produce transverse forces of similar magnitude. Figure 
5 shows that this is also true for plates having a thickness in the 
neighborhood of 41. There is therefore no preferred thickness with re- 
spect to favoring applied field gradients over temperature or com- 
position gradients. 

In temperature dependent materials in which the domain tends to 
move towards the high temperature direction the domain will tend to 
follow a laser beam initially placed at its center. In materials having 
the opposite temperature characteristic, the domains may be pushed 
by a laser beam. Gradients in thickness or composition may be used 
to define domain tracks in order to increase margins with respect to 
spurious field or temperature gradients. An immediately useful ap- 
plication of equation (18) is in the calculation of the effect on device 
margins of the heat produced by domain generators and detectors. 
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APPENDIX 

The Internal Magnetostatic Energy Function 

Various representations of the internal magnetostatic energy func- 
tion, I(d/h), are given in this appendix. Since F{d/h) is non-negative 
(Ref. 3, Fig. 3 and equation 138), the integral definition of I(d/h) (2) 
implies that I(d/h) is positive and monotonic increasing as can be 
seen in Fig. 1. (Neither function is defined in the present case for 
negative values of the argument.) The definition of I(d/h) also 
implies the differential equation, dl(x)/dx = F(x), with the boundary 
condition 1(0) = 0. The internal magnetostatic function may be 
written as 



±d r 

3irhl 
where as in equations (84) and (85) of Ref. 3 
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and 



V(x) - 2(x + l)- i x( r ^) 



(19) 



(20a) 



(20b) 



where K(m) and E(m) are the complete elliptic integrals of the first 
and second kind respectively and the argument is in the m of Ref. 8. 
Differentiation of equation (19) with respect to d/h and combining 
terms using equations (86) and (87a) of Ref. 3 verifies that the 
derivative of I (d/h) is indeed F(d/h) (equation 138 of Ref. 3). Sub- 
stitution of the series expansions of U and V (equations 96 and 97 of 
Ref. 3) verifies that 7(0) = 0. 

Writing I (d/h) in terms of the Li function (Ref. 3 Sec. A.4) elimi- 
nates the negative power of d/h, 
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(21) 



and the expression in terms of F(d/h) and S (t (d/h) (equations 138 
and 139 of Ref. 3) is also sometimes useful, 
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The power series expansions of I(d/h) may be determined either 
by substituting the power series expansions for U, V and L x (Ref. 3, 
Sections A.3 and A.4) into equation (21) or by term-by-term integra- 
tion of F(d/h) with the integration constants being determined from 
the lowest order terms of the power series expansions of U and V. The 
result is 
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